Weighted Density Approximation Description of Insulating YH3 and LaHs 
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, Density functional calculations within the weighted density approximation (WDA) are presented 

I ■ for YH3 and LaHa. We investigate some commonly used pair-distribution functions G. These cal- 

culations show that within a consistent density functional framework a substantial insulating gap 
I can be obtained while at the same time retaining structural properties in accord with experimen- 

^ . tal data. Our WDA band structures agree with those of GW approximation very well, but the 

' calculated band gaps are still 1.0-2.0 eV smaller than experimental findings. 

CNi : 

, The rare earth trihydrides, MH3, where M = Y, La - Lu, are the subject of considerable practical and theoretical 
= interest. The rare earths absorb hydrogen readily and in the trihydride form have among the largest hydrogen to metal 
[ ratios of the elemental hydrides The fully hydrogenated materials are insulating, while those with less hydrogen 
I are metallic, and in thin film form they can be used as switchable mirrors It is important to demonstrate a correct 
description of these trihydrides in first-principles approaches in order to validate calculations for more complex metal 
hydrogen systems. However the standard first-principles calculations based on density functional theory (DFT) within 
the local density (LDA) and generalized gradient (GGA) approximations incorrectly predict metallic behavior for the 
, rare earth trihydrides. Here we explore application of the weighted density approximation (WDA). We find a WDA 
^ ' that works for these materials (and others), indicating that the gap problem here is not due to unusual correlations 
Q I or quasiparticle corrections, but is simply a problem with the LDA and GGA exchange correlation (xc) functionals. 
O The rare earth dihydrides, MH2, are metallic. The trihydrides on the other hand are insulating with optical gaps 

^ ^ of the order of 2 eV, with a metal insulator transition at a hydrogen content correlated with rare earth ionic radius 
• |3- This insulating gap of the trihydrides may be understood in terms of simple picture, in which H forms a negative 
^\ ' ion, H~, fully ionizing the metal atoms, thus forming a gap between H bands and higher lying metal sd conduction 
bands. This simple picture is supported by infrared and photoemission measurements 3, 4j. It is consistent with 
the fact that all the rare earth trihydrides are insulators despite differences in the details of their crystal structures. 
These range from the fee derived structure of LaHs to the not fully solved hexagonal structure of YH3 0, 0, 0, • 

Computations within the LDA predict a metallic band structure with substantial overlap between H and metal 
derived bands jSi, iSi, ilfll lUl . This severe gap problem could be due to several possible causes: (1) use of an incorrect 
crystal structure; (2) an extreme case of LDA error in predicting gaps; (3) unusual correlation effects; or (4) self- 
interaction errors in LDA. Early attempts at resolving this discrepancy focused on crystal structure. Indeed it was 
shown that with optimized crystal structures, DFT calculations could yield insulating band structures for YH3 15,1^. 
Unfortunately, these relaxed structures are not consistent with diffraction measurements, and such an explanation 
would not explain the similar insulating properties of the other AfHa compounds. 

DFT is a theory of ground state properties, and the eigenvalues of the Kohn-Sham equations cannot be equated 
to quasiparticle excitations 1A\ ■ Nevertheless, any non-pathological accurate approximation to the exact DFT must 
predict an insulating Kohn-Sham band structure, because static responses and other ground state properties depend 
on the presence or absence of a gap. For example, LO-TO splitting would be present in the phonon spectrum of ionic 
insulator like YH3, but not in a metal. 

Recent studies on rare earth trihydrides have considered electron correlation effects 0, 0, 0, 0, 0| , particularly 
electron correlations on the H sites. However, since there are no partially occupied d 01 f orbitals in LaHs and YH3, 
nor is there evidence that these materials are near a quantum critical point, these hydrides would represent a rather 
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novel class of strong correlated materials. Calculations done within the GW approximation show that band gaps that 
are reconcilable with experiment are obtained, both with the cubic BiF3 structure (LaHs) and the hexagonal H0D3 
structure (a simplified YH3 structure) 0, IS ■ The GW method approaches the problem as being most simply 
viewed as due to the treatment of excited states, but if the problem is due to self-interaction, it is fundamentally 
a ground state problem. The LDA makes large errors in description of the one-electron hydrogen atom because of 
incomplete cancellation of the Hartree and xc self-interactions. One way to remove the spurious self interactions in 
LDA and GGA calculations is the so-called self interaction correction 20]. In this approach an orbital dependent 
self interaction is calculated at each step and used to create an orbital dependent single particle Hamilitonian. This 
method suffers from basis set dependence, e.g. the computed correction is zero for a Bloch state representation. We 
explore whether one can remain within conventional ground state DFT and obtain the proper electronic structure of 
rare earth trihydrides. 

The weighted density approximation (WDA) |21|, |23| has the exact non-local form of DFT expression of xc energy, 

E^,[n] ^ j j y"y [g.e(r,rO - l]dvdv' . (1) 

Since the exact pair-distribution function .gxc(r, r') is elusive, a model WDA function G is made, 

g,e(r,r')-l = G[|r-r'|,n(r)], (2) 
where the weighted density n can be determined from the xc hole condition: 

j n(r')G[|r - r'|, n(r)]dr' = -1. (3) 

This assures no self-interaction in a one-electron system. The WDA also gives the correct result for the uniform electron 
gas because it reduces to the LDA in this case. A seemingly natural choice of G would be of a uniform system, but 
Mazin and Singh [23| emphasized, based on the xc kernel, that G of the uniform electron gas is not optimal. People 
advanced several model forms of G, such as the Gunnarsson- Jones (G-J) ansatz, G^^''(r, n) = c{l — exp(— [j]^'')}, 
with fc = 5; and the Gritsenko et al. (GRBA) ^ ansatz, G^^^^(r, n) = c • exp(-[J]''), with k = 1.5. In this letter, 
we also used a homogeneous type of G and another G-J type for = 4 (G-J4), which was tested by Rushton et al. 
for silicon 0| . To avoid the unphysical WDA inter-shell interaction, we used the shell-partition method 0, |2^ . 

As mentioned before, although the exact structure of hexagonal YII3 is still under investigation, its insulating 
property is believed to be an electronic property, i.e., not sensitive to the crystal structure (301 jsij . So we studied 
a simplified hexagonal LaF3 structure with a Y2H6 unit cell |llj |. in which Y atoms form a hexagonal close-packed 
(hep) structure. The measured lattice constants are a — 3.672 A and c/a — 1.81. We optimized the H positions using 
LDA and WDA (with four types of G), and all the results are consistent with the previous LDA result ^11]. We also 
optimized the equilibrium volume with the calculated atomic positions and fixed c/a ratio. Table shows that the 
WDA volumes are only about 2-3% deviated from experimental data. 

To determine the band structure of the hexagonal YII3, We used the experimental volume and al2xl2x8 special 
k-point mesh. The calculated bands of the LDA and WDA (with the G-J4 type of G) are shown in Fig. ^ The LDA 
predicts a direct band gap of 1.0 eV at F and an overlap of 0.93 eV between F and K, which leads to a semi-metal. 
Our present LDA results are a little different from the previous (0.6 eV for the F gap and 1.3 eV for the overlap 
[ill) because we used the Hedin-Lundqvist xc function, whereas they used the Wigner ansatz. On the other hand, 
the WDA predicts a direct F gap of 2.2 eV, and opens up a fundamental gap of 0.41 eV between F and K. Thus the 
WDA predicts a substantial insulating gap. Our WDA band gaps are pretty close to the GW calculations 17] of 0.6 
eV for the fundamental gap and 2.9 eV for the direct F gap. In the Brillouin zone of the more complicated tripled 
H0D3 structure of YH3, this indirect F-K gap will be folded to form a direct gap at F point, and they are roughly of 
the same order of magnitude [iJi E| . 

We also calculated the WDA bands with the other three types of G, and the results are shown in Table ITTl We found 
the WDA with all these G functions gives better band gaps than the LDA. Except for the G-J4 form, the GRBA 
form is superior to the other two, because it predicts a smaller band overlap between F and K. These results illustrate 
that the WDA band structure is sensitive to the choice of G. We compare these pair-distribution functions by drawing 
them together in Fig. [3 where n = 3/{4TTr^), and in the WDA calculations of YH3, 0.8 < < 2.3. We found the 
G-J4 form has the shortest range, and the range of GRBA ansatz is shorter than the other two. Our calculations 
indicate that the WDA with shorter range G gives larger band gaps. Interestingly, such functions, particularly GRBA, 
yield static responses closer to Monte Carlo results for the uniform electron gas |24| . 

We performed LDA and WDA (with G- J4 only) calculations on cubic B1F3 structure of YH3 and Lall3 , in which the 
Y or La atoms constitute a face-centered cubic (fee) lattice and H atoms are located at the tetrahedral and octahedral 
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sites. Recently the insulating fee YHs was stabilized by MgH2 |3l|. There is a large discrepancy of the measured 
volumes of fee YH3. Ahuja et al. |0 expected 32.2 A'^, while van der Molen et al. [sj extrapolated 36.9 A'^ and 
Gogh et al. I^D] obtained 38.0 A^. Our LDA and WDA equilibrium volumes are 32.4 A^ and 35.2 A'^ respectively. 
We used the experimental volume of 36.9 A^. For LaHs the measured volume 44.2 A^ is close to LDA and WDA 
equilibrium volumes. We used a 12 x 12 x 12 special k-point mesh. 

In Fig. 121 are the calculated band structures of cubic YH3 and LaHa. The LDA predicts that both of them semi- 
metals with direct overlaps of 1.0 eV for YH3 and 0.46 eV for LaHa at F point, which agree with previous calculations 
[Tglls^ . Contrary to the LDA, the WDA predicts that both trihydrides are semiconductors with direct F gaps of 0.19 
eV and 0.80 eV for YH3 and LaHa respectively. However the fundamental gaps are about 0.1 eV smaller, because 
the lowest conduction bands at L is slightly lower than that at F. Our WDA band gaps are in good agreement with 
previous GW calculations. Alford et al. |32| predicted a F gap of 0.2-0.3 eV for YH3 and 0.8-0.9 eV for LaHa, and 
Chang et al. obtained 0.8 eV for LaHa. We found that the hexagonal and cubic YH3 have similar small gaps, 
suggesting that the insulating properties are of electronic origin rather than structural. 

Another interesting feature of the WDA band structures is that compared with the LDA, the overall valence 
bandwidth hardly changes. This is also consistent with previous GW calculations [TtL IT^ l32| , but in contrast to 
the many-body model Hamiltonian calculations which show a large valence bandwidth decrease |15|. It confirms the 
previous conclusion [TtI Is^ that rear-earth trihydrides are just band insulators, and that strong correlations are not 
required. 

Although our DFT-WDA band gaps of YII3 and LaHa agree with previous GW quasiparticle calculations very well, 
the measured optical gap is ~2.8 eV for YH3 and ~1.8 eV for LaHs, which are still 1-2 eV larger. Gelderen et al. 
[f^ argued that the fundamental band gap of hexagonal YH3 is not the measured optical gap due to the forbidden 
transition between the highest valence and the lowest conduction bands, and they concluded that the measured 
value corresponds to the direct gap at F of the hexagonal LaFs structure. However the real hexagonal YH3 is more 
complex and it may have low symmetry so that the dipole matrix element is not zero. Our DFT-WDA method may 
underestimate band gaps due to the DFT gap discontinuity at the fermi level j3^ . 

In conclusion, we have studied the band structures of the hexagonal YH3 and cubic YH3 and LaH3 using the first- 
principles WDA method. The calculations showed that YH3 and LaH3 are both semiconductors, and this insulating 
property does not rely on the detailed crystal structure. The calculated band gaps agree with previous GW calculations 
very well, and they are also consistent with experiment. Our WDA calculations also predicted good ground-state 
structures. The success of WDA may be attributed to the absence of self-interaction in H Is sates. 
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TABLE I: With the calculated atomic positions and fixed c/a ratio, the optimized equilibrium volumes (in A'') of hexagonal 
YH3 (Y2H6). Here a, b, c and d denote uniform, GRBA, G-J and G-J4 forms of G respectively. Numbers in parentheses are 
the percentage deviations from experiment. 
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TABLE IL Calculated LDA and WDA (a, b, c and d as seen in Table|IJ fundamental and V gaps (in eV), compared with GW 
results 17]. A negative band gap means band overlap. 
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FIG. 2: Four types of the pair-distribution function G(r, n) with Ts — 0.8 and — 2.3. 
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FIG. 1: LDA and WDA (with the G-J4 type of G) electron bands of the hexagonal YII3. 
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FIG. 3: LDA and WDA (with the G-J4 type of G) electron bands of the cubic YH3 and LaHs. 



